Stable propagation of a modulated particle beam in a crystal channel 
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The propagation of a modulated beam of charged particles in a planar crystal channel is investi- 
gated. It is demonstrated that the beam preserves its modulation at sufficiently large penetration 
depths to ensure the feasibility of using a crystalline undulator as a coherent source of hard X rays. 
This finding is a crucial milestone in developing a new type of lasers radiating in the hard X ray 
and gamma ray range. 
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In this Letter we study for the first time the evolution 
of a modulated particle beam in a planar crystal channel 
and demonstrate that it preserves its modulation at suf- 
ficiently large penetration depths to ensure the feasibility 
of using a crystalline undulator as a coherent source of 
hard X rays. Solving this problem is of crucial impor- 
tance in the theory of the Crystal Undulator based Laser 
(CUL) [1, 2, 3] — a new electromagnetic radiation source 
in hard X and gamma ray range. 

Channeling takes place if charged particles enter a sin- 
gle crystal at small angle with respect to crystallographic 
planes or axes [4]. The particles get confined by the in- 
terplanar or axial potential and follow the shape of the 
corresponding planes and axes. 

A single crystal with periodically bent crystallographic 
planes can force channeling particles to move along nearly 
sinusoidal trajectories and radiate in hard X and gamma 
ray frequency range (see Fig. 1). The feasibility of such a 
device, known as the 'crystalline undulator', was demon- 
strated theoretically a decade ago [i] (further develop- 
ments as well as historical references are reviewed in ["i]). 
Its experimental study is on the way within the PECU 
project [()]. 




Periodically bent 
^crystallographic 
\plane 



FIG. 1: Schematic representation of the crystalline undulator. 

The advantage of the crystalline undulator is in ex- 
tremely strong electrostatic fields inside a crystal which 
are able to steer the particles much more effectively than 
even the most advanced superconductive magnets. Due 



to this fact, such an undulator can radiate powerful elec- 
tromagnetic waves in the hard X and soft gamma ray 
range, where conventional sources with comparable in- 
tensity are unavailable [7]. 

Even more powerful and coherent radiation will be 
emitted if the probability density of the particles in the 
beam is modulated in the longitudinal direction with the 
period A, equal to the wavelength of the emitted radia- 
tion. In this case, the electromagnetic waves emitted in 
the forward direction by different particles have approx- 
imately the same phase [8]. Therefore, the intensity of 
the radiation becomes proportional to the beam density 
squared (in contrast to the linear proportionality for an 
unmodulated beam). This increases the photon flux by 
orders of magnitude relative to the radiation of unmodu- 
lated beam of the same density. The radiation of a mod- 
ulated beam in an undulator is a keystone of the physics 
of free-electron lasers [!)]. It can be considered as a clas- 
sical counterpart of the stimulated emission in quantum 
physics. Therefore, if similar phenomenon takes place in 
a crystalline undulator, it can be referred to as the lasing 
regime of the crystalline undulator. 

The feasibility of CUL radiating in hard X ray and 
gamma ray range was considered for the fist time in [1, 2]. 
Recently, a two-crystal scheme, the gamma klystron, has 
been proposed [•')]. 

A simplified model used in the cited papers assumed 
that all particle trajectories follow exactly the shape of 
the bent channel. In reality, however, the particle moving 
along the channel also oscillates in the transverse direc- 
tion with respect to the channel axis (see the shape of 
the trajectory in Fig. 1). Different particles have dif- 
ferent amplitudes of the oscillations inside the channel 
(Fig. 2, upper panel). Similarly, the directions of parti- 
cle momenta in (xz) plane are slightly different (Fig. 2, 
lower panel). Even if the speed of the particles along 
their trajectories is the same, the particles oscillating 
with different amplitudes or the particles with different 
trajectory slopes with respect to z axis have slightly dif- 
ferent z components of their velocities. As a result, the 
beam gets demodulated. An additional contribution to 
the beam demodulation comes from incoherent collisions 
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FIG. 2: Due to different amplitudes of channeling oscillation 
(upper panel) and different momentum directions in the (xz) 
plane (lower panel), the initially modulated beam gets de- 
modulated. Open and filled circles denotes the same particles 
at the crystal entrance and after traveling some distance in 
the crystal channel, respectively. 



of the channeling particles with the crystal constituents. 

In the case of an unmodulated beam, the length of 
the crystalline undulator and, consequently, the maxi- 
mum accessible intensity of the radiation is limited by 
the dechanneling process. The channeling particle grad- 
ually gains the energy of transverse oscillation due to 
collisions with crystal constituents. At some point this 
energy exceeds the maximum value of the intcrplanar po- 
tential and the particle leaves the channel. The average 
penetration length at which this happens is known as the 
dechanneling length. The dechanneled particle docs not 
follow the sinusoidal shape of the channel any more and, 
therefore, docs not contribute to the undulator radiation. 
Hence, the reasonable length of the crystalline undulator 
is limited to a few dechanneling lengths. A longer crys- 
tal would attenuate rather then produce the radiation. 
Since the intensity of the undulator radiation is propor- 
tional to the undulator length squared, the dechanneling 
length is the main restricting factor that has to be taken 
into account when the radiation output is calculated. 

In contrast, not only the shape of the trajectory but 
also the particles positions with respect to each other 
along z axis are important for the lasing regime. If this 
positions become random because of the beam demod- 
ulation, the intensity of the radiation drops even if the 
particles are still in the channeling mode. Hence, it is 
the beam demodulation rather than dechanneling that 
restricts the intensity of the radiation of CUL. Under- 
standing this process and estimating the characteristic 
length at which this phenomenon takes place is, there- 
fore, a cornerstone of the theory of this new radiation 
source. 

Let us consider the distribution /(t, z; ^, Ey) of the 



beam particles with respect to the angle between the 
particle trajectory and axis z in the (xz) plane ^ = 
avcsmpx/p ~ Px/p and the energy of the channeling os- 
cillation Ey = Py/2E + U{y) [11]. Here p, p^ and py 
are, respectively, the particle momentum and its x and y 
components, U (y) is the interplanar potential, and E is 
the particle energy (we will consider only ultrarelativis- 
tic particles, therefore E k. p). It can be shown that 
the evolution of this distribution in the crystal channel 
with the time t and the longitudinal coordinate (pene- 
tration depth into the crystal) z can be described by the 
following differential equation of Fokker-Planck type: 
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Here Dq is the diffusion coefficient that is dominated by 
the scattering of the beam particles by lattice electrons. 
The particle longitudinal velocity averaged over the un- 
dulator period, Vz, is given by 
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If the beam is periodically modulated (bunched) the 
distribution can be represented as a Fourier series: 
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with g*{z;£,,Ey) = g-j{z; ^, Ey) to ensure the real value 
of the particle distribution. Since Eq. (1) is linear, it 
is sufficient to consider only one harmonic. Substituting 
fit, z; ^, Ey) = g{z] ^, Ey) exp^iut) one obtains 
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To simplify this equation, we make the substitution 
g{z;^,Ey) = exjp {—iujz) g{z; ^, Ey) and assume that the 
variation of g{z;^,Ey) within the modulation period is 
small: dg/dz ^ ujg{z; Ey). This allows us to neglect 
the terms (1 — Vz)dg/dz while keeping the terms (1 — 
Vz)ujg{z; ^, Ey). The resultant partial differential equa- 
tion for g{z; Ey) can be solved by the method of sepa- 
ration of variables. Putting g{z; Ey) = Z{z)'E.{£^)£{Ey), 
we obtain a set of ordinary differential equations: 
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where Cz , Q and Cy do not depend on any of the variables 
z, ^ and Ey and satisfy the condition 



Cz = Ce 



(8) 
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Eq. (5) has the form of the Schrodinger equation for 
the harmonic oscillator. Its eigenvalues and eigenfunc- 
tions are, respectively, 



Ci,n = -{l + i)\l^(n+l] , n = 0,l,2,. 



E 
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and 



Here . . ) arc Hermitc polynomials. 

Eq. (6) can be reduced to the Laguerre differential 
equation, so that its solution can be represented as 



(11) 

where . . ) is the Laguerre function[i2] and Vk is re- 
lated to the eigenvalue Cy^k via 



Cy,k = -^-^\l^{2iyk + l), fc = 1,2,3,... (12) 

The eigenvalues can be found by imposing the bound- 
ary conditions: the density of the channeling particles 
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becomes zero if the energy of channeling oscillations Ey 
equals to the interplanar barrier Umax- 
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Equation (13) has to be solved for Vk- (the subscript k 
enumerates different roots of the equation) and the result 
has to be substituted into (12). 

It is convenient to represent the eigenvalues in the form 



(14) 



Here Ld = 4J7inax/(jo,i^o) is the dechannchng length 
[10] (io,fc is fc-th zero of the Bessel function Jo(e)), 
6'l = ^j2Ui-ae,x/ E is Lindhard's angle. We introduced 
the parameter 



(15) 



where A = 27r/w is the spatial period of the modulation. 
The functions ak{n) and Pkii^) are to be found by solving 
numerically Eq. (13) combined with (12). 
Using (8), one finds the solution of Eq. (7): 
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Hence, the solution of Eq. (4) is represented as 
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where the coefficients d^fc are found from the particle dis- 
tribution at the entrance of the crystal channel. Due to 
the exponential decrease of Zn,k{z) with z (see (16)), the 
asymptotic behaviour of g{z;S,,Ey) at large z is domi- 
nated by the term with n = and fc = 1 having the 
smallest value of the factor [q;/c(k) 4- (2n + 1)\/k/jo,i] in 
the exponential. Therefore, at sufficiently large penetra- 
tion depths, the particle distribution depends on z as 
g{z;^,Ey) cx cxp (— z/idm ~ i'^-'/uz -z) where Ldm is the 
newly introduced parameter — the demodulation length: 
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and Uz is the phase velocity of the modulated beam along 



the crystal channel 
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This parameter is important for establishing the reso- 
nance conditions between the undulator parameters and 
the radiation wavelength. It will be analyzed elsewhere. 

In this Letter we concentrate our attention on the de- 
modulation length. This parameter represents the char- 
acteristic scale of the penetration depth at which a beam 
of channeling particles gets demodulated. Fig. 3 presents 
the dependence of the ratio Ldm/id on the parameter 
K. It is seen that the demodulation length approaches 
the dechanneling length at k < 1. This is an impor- 
tant result, since it means that the demodulation pro- 
cess docs not put additional restrictions on the undula- 
tor length and, therefore, a strong lasing effect can be 
expected in the crystalline undulator fed by a modulated 
beam. On the contrary, the ratio noticeably drops for 
K > 10. Therefore, further investigations are needed to 
clarify whether the lasing regime is feasible for these k. 



4 









- 




- 




Motion in x-direction only \ '\ 






Motion in y-direction only ^s/'**-' 






Both contributions are included 





0.01 0.1 1 10 100 1000 



K 

FIG. 3: The ratio of the demodulation length Ldm (18) to the 
dechanneling length Ld versus the parameter k (15). See text 
for details. 

The dependence of the parameter k on the energy of 
the emitted photons, fku ~ 2T:h/X, is shown in Fig. 4. 
The calculation was done for 1 GcV[ ' '] positrons using 
the formula for the dechanneling length from [•^), 10]. 
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FIG. 4: The parameter k (15) versus the photon energy hui 
for different crystals and crystallographic planes. 

As one sees from the figure, there exist such crystal 
channels that k — 1 corresponds to huj ~ 100 — 300 keV. 
Therefore, CUL is most suitable for the application in 
this photon energy range. Using crystalline undulator for 
the emission of softer photons huj < 10 keV is prevented 
by very strong absorption of the radiation in the crystal, 
while the upper limit on the photon energy is set by the 
decreasing value of Ldm ■ 

It is instructive to study the influence of the particle 
motion in x and y direction on the demodulation length 
separately. Replacing aii^n) in (18) with unity means 
neglecting the motion in the y direction, while omitting 
the second term in the denominator ignores the motion 
in X direction. One sees from Fig. 3 that it is mostly the 
motion in x direction that suppresses the demodulation 
length at k < 10, while the influence of channeling os- 
cillations is negligible. This suggests the idea that axial 



channeling, i.e. when motion in both x and y directions 
has the nature of channeling oscillations, might be more 
suitable for the lasing regime in the range hoj ^ 1 MeV. 

Iir conclusion, we have studied the propagation of a 
modulated particle beam in a planar crystal channel. It 
has been demonstrated that the beam preserves its mod- 
ulation at the penetration depths which are sufficient for 
using the crytalline undulator as a source of coherent ra- 
diation with the photon energy of hundreds keV. 

We restricted our analysis to the beam demodulation 
in a straight crystal. The influence of the periodic bend- 
ing on the demodulation length has to be studied at the 
next step. Another important milestone in the theory 
of CUL would be developing suitable methods of beam 
modulation. 
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